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Abstract 

Active thermography is a non-destructive testing technique to detect the internal structure of a heat 
conductor, which is widely applied in industrial engineering. In this paper, we consider the problem of 
identifying an unknown cavity with Robin boundary condition inside a heat conductor from boundary 
measurements. To set up the inverse problem mathematically, we hrst state the corresponding forward 
problem and show its well-posedness in an anisotropic Sobolev space by the integral equation method. 
Then, taking the Neumann-to-Dirichlet map as mathematically idealized measured data for the active 
thermography, we present a linear sampling method for reconstructing the unknown Robin-type cavity 
and give its mathematical justihcation by using the layer potential argument. In addition, we analyze the 
indicator function used in this method and show its pointwise asymptotic behavior by investigating the 
reflected solution of the fundamental solution. From our asymptotic analysis, we can establish a pointwise 
reconstruction scheme for the boundary of the cavity, and can also know the distance to the unknown 
cavity as we probe it from its inside. 

Keywords. Inverse boundary value problem; Heat equation; Robin boundary condition; Reconstruc¬ 
tion scheme; Asymptotic behavior. 
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1 Introduction 


Active thermography is a widely used non-destructive testin g te chnique in industrial engineering, which aims 
to detect the internal information of a heat conductor [?, 21, 36, 371. The principle of active thermography is 


as follows. If there is an anomaly inside the conductor, it affects the propagation of the heat flow inside the 
conductor and as a result it also directly affects the temporal behavior of the surface temperature distribution. 
By measuring the surface temperature distribution and by solving some inverse problem, the information 
on the anomaly can be calculated. The information we want to know are the size, location and shape of 
anomalies and their physical properties such as heat conductivities. The measurement of active thermography 
is a non-contact, very fast and large area measurement which is conducted by injecting a heat flux to the 
conductor by a flash lamp or heater and measuring the corresponding distribution of temperature on the 
surface of the conductor by an infrared light camera. 
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In this paper, we want to recover an unknown cavity with Robin boundary condition inside a heat 
conductor via active thermography. To begin with, we give the mathematical formulation as follows. Let 
fl C K" (n = 2, 3) be a heat conductor and U a cavity with Robin boundary condition embedded in fl. We 
assume that the boundaries dfl and dU of Q and U, respectively, are of class C^. Injecting a heat flux / on 
dfl over some time interval (0, T), the corresponding temperature distribution u(x, t) in fl \ D x (0, T) can 
be modeled by the following initial-boundary value problem: 

{ (dt-A)u = 0 in (n\B) X (0, T), 
dvU — Am = 0 on dD x (0, T), 
d^u = f on dQ X (0, T), 

M = 0 at t = 0, 


where A = A(a;) € {9D) is the real-valued impedance and i/ on dD (or dfl) is the unit normal vector directed 

into the exterior of D (or D). We will show that the initial-boundary value problem (ll.l|l is well-posed in a 
suitable Sobolev space. Then we idealize a set of many pairs of the heat flux / and the corresponding surface 
temperature distribution u\gux{o,T) which is the measured data for active thermography as the Neumann- 
to-Dirichlet map Ajj given by Ajj : / i—>■ u\gux{o,T)- Thus, our inverse problem for thermography is to 
reconstruct D from Ao- 

When D is a usual cavity with Neumann boundary condition or an inclusion, the corresponding inverse 
problem has been extensively studied. In [s, 14-16|, the uniqueness and stability estimate are established. 
In , Newton-type iteration algorithms based on domain derivatives are studied. As for non-iterative 

reconstruction schemes, we can consult the papers 12, 13, 22, 23, 25, [i^ 32, 34, and the references 


therein, where the reconstruction schemes called the dynamical probe method and the enclosure method 
are extensively studied. A reconstruction scheme for unknown cavities via Feynma-Kac type formula is also 
proposed in [ 23 . Recently, the authors established a linear sampling-type method for the heat equation to 
identify unknown cavities with Neumann boundary condition |^ . This method was extended to the inclusion 
case in 


35|. 



In this paper, we are concerned with the reconstruction of an unknown cavity where a Robin boundary 
condition is prescribed. Some results on uniqueness and stability for this inverse problem can be found in 
0,0. Assuming that D is given, the determination of the Robin coefficient A from measured data was 
also considered; see 1^ 271 and the references therein. In this work, assuming that both the geometrical 
information of the cavity D and the Robin coefficient A on its boundary are unknown, we establish a linear 
sampling method to identify the unknown Robin-type cavity D from the boundary measurements An . As 
we know, the linear sampling method was originally proposed for inverse scattering problems in 

was further investigated from both the theoretical and numerical aspects; for example, see 0 , 0 __ 

31 1 Hi 0 0 • Roughly speaking, this method in the heat equation case is based on the characterization of the 
approximate solvability of the so-called Neumann-to-Dirichlet map gap equation (A^j —Agj^* = r°(a:, t; y, s), 
where A 0 is the Neumann-to-Dirichlet map when there are not any cavities inside D, and r''(a:, t; y, s) is the 
Green function for the heat operator dt — A in Dx (0, T) with homogeneous Neumann boundary condition on 
its boundary. By giving this characterization, we define a mathematical testing machine called an indicator 
function to reconstruct the boundary of D. When we probe it from the inside of D, we don’t need to let the 
discrepancy of the Neumann-to-Dirichlet map gap equation tend to zero. 

The new ingredients of this paper consists of two parts. First, we prove the well-posedness of the forward 
problem (HH) in an anisotropic Sobolev space by the integral equation method. This is necessary for our 
purpose to investigate the inverse problem using layer potential argument, although the well-posedness can 
also be justified in the usual function space W(0, T) by the argument in [^. Second, compared with our 


previous works in 


H 


for the Neumann cavity and inclusion cases, we provide a further investigation of 


the linear sampling method for the heat equation and show the asymptotic behavior of the indicator function 
by carefully analyzing that of the reflected solution of the fundamental solution. Since the reflected solution is 
the compensating term of the Green function for the related initial-boundary value problem, we actually give 


2 






























a pointwise short time asymptotic behavior of the Green function near dD. As a consequence, by observing 
the asymptotic behavior of the indicator function, we can have the followings: 

(i) pointwise reconstruction of the shape and location of an unknown cavity; 

(ii) information about the distance to the boundary of the cavity. 

The more precise meaning of the above item (ii) is that we can know the distance to dD when we probe it 
from inside D by using the indicator function. 

The rest of the paper is organized as follows. In Section [2l we show the well-posedness of the forward 
problem in an anisotropic Sobolev space by the integral equation method. Then, in Section [3l we present 
a mathematical justification of the linear sampling method, while the asymptotic behavior of the indicator 
function is provided in Section 31 Finally, in Section [Sj we give some concluding remarks. 


2 The integral equation method for the forward problem 

In this section, we show the unique solvability of the forward problem by the integral equation method. Then 
we explicitly define the Neumann-to-Dirichlet map A^i and state our inverse problem. Since our analysis for 
the inverse problem is based on the layer potential argument, it is necessary to establish the well-posedness 
of the forward problem in an anisotropic Sobolev space, instead of the usual function space VF(0, T). 

Let us start by introducing the anisotropic Sobolev spaces. For p, q > 0 we define 

X K) := HP{W^)) n L2(K")). 

For p, q < 0 we define the space by duality x R) := x K))^ Throughout this paper, 

we denote X x (0, T) and dX x (0, T) by Xt and {dX)T, respectively, where AT is a bounded domain in K” 
and dX denotes its boundary. By HP’‘^{Xt) we denote the space of restrictions of elements of x R) 

to Xt- The space H^'‘^{{dX)T) is defined analogously. We also introduce the following function spaces: 

{Xt) := G {X x (—oo, r))| u{x, t) = 0 for t < o|, 

H^’^Xt; dt - A) := {u e H^'^{Xt)\ {dt - X)u G L^{Xt)} - 

Then our forward problem is formulated as follows. 

For-ward problem: Given f € H 2’ I ((dlljr) and A G C^{dD), find a unique solution u G ((fl \ 
D)t) to the problem (11.11) . 


Denote by 


r(a;, t; y, s) := 


_i_exp ( 

0 , 


t > s, 
t < s 


the fundamental solution of the heat operator dt — A. For convenience, we sometimes write it as F(y j;)(x, t). 
Define the following heat layer potentials: 

{V,j(p){x, t) ■-= [ [ r(x, t; y, s)(p{y, s)da{y)ds, {x, t) G Sj x (0, T), 

Jo JSi 

{Kijip){x, t) ■= J s)da{y)ds, {x, t) G Sj x (0, T), 

{Nijip){x, t) := [ ( I’ ^\ {y, s)da{y)ds, (x, t) G Sj x (0, T), 

Jo JSi av{x) 

{W,jip){x, t) ■-= -tA-t f [ ^^ <p{y, s)da{y)ds, {x, t) G Sj x (0, T). 

onx) Jo Jsi du{y) 
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In this paper we take i, j = 1, 2 with = dD and S 2 = dfl. 

We now show the unique solvability of (HI) in ((n \ D)t) by the integral equation method. 

Theorem 2.1 Suppose that A £ C^{dD) and 0 < Ai < A < A 2 , where Ai and A 2 are two given constants. 
Then, there exists a unique solution u £ ((O \ D)t) to the initial-boundary value problem (11.11) for any 

f€H-h-i{{dn)T). 

Proof. To prove the uniqueness of solutions in \ D)t), we first establish the uniqueness in the 

function space W (0, T) defined by 

W(0, T):=^w: we L^((0, T); \ D)), dtw £ L" ((0, T); \ D))') } . 

Consider the bilinear form associated with (H 

b{u, V, t) ■.= / Vu-\7vdx-\- / Xuv da{x) 

Jn\D Jan 

for u, V e D) and t £ [0, T], It is easy to verify that 

\b{u, V, t)| < '"ll^llffi(n\'D) ll^ll//i(n\'D)5 

where c is a positive constant. Then the uniqueness of solutions to (H in tT(0, T) can be justified in a 
standard way; see, for example, [^, §26]. Since any solution u £ ((fl \ D)t) to {dt — A)u = 0 with zero 
initial condition is in IT(0, T), we also have the uniqueness in i /^’2 ((fl \ D)t). 

To prove the existence of the solution in \ D)t), we make use of the integral equation method 

ll| . By Green’s representation theorem for (11.111 . we have 


L{x,t)= f f _ r(a:, t; y, da{y)ds - 

Jo Jd{n\D) ov[y) 

= -/7 nx,t;y,s)^-^daiy)ds + r [ 
Jo JdD ov{y) Jq Jq 

+ /7 n..t,y.s)?^My)'ls- f j 

Jo Jan ov{y) Jq Jq 


dT{x, t; y, s) 


0 Jd{Q\D) 9 v{y) 
* '■ ar(a;, t; y, s) 


u{y, s)da{y)ds 


dv{y) 

dT{x, t] y, s) 
0 Jan 9 i^{y) 


u{y, s)da{y)ds 


u{y, s)da{y)ds. (2.1) 


T / \ ^ 

Set (ui, U2, M3) = {duu\(^gD)TJ 'J'\iaD)T^ '^l(ao)T) • Then, using the jump relations of layer potentials, we 

deduce the following system of boundary integral equations: 


K 2 


-i/ + 7 Vii Wii+XI -IT21 


Ml 

M 2 


Tl2 —Ki 2 hi + K22 J V M3 


V..f \ 

= I N2if . 
1 ^ 22 / / 


( 2 . 2 ) 


Define the space % .= H 2’ i {{dD)T) x H^’l [[dD)T) x (((?D)t) with the norm 


||(mi,M2,M3)7|«:= IKII; . + +||m3|H ^ 

' H 2' 4((0D)t) H2'A[[dD)T) H 2'A [[dVtjT) 


1/2 


for (mi, M 2 , £ H. We show the unique solvability of the system ( 12 . 21 ) in %. To accomplish this, we 

define the operator 

Vn \I - Kn ' 


A:= 


-XI + Nil Wii + XI 


( 2 . 3 ) 
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and the space Ho ■= H 2 ’ 4 ((i9D)t) x H'^'i{{dD)T) equipped with the norm 

1/2 


||(ui, M2)' liwo := lluir 1 1 +||M2|rii 


for (mi, U2Y £ Ho- 


Let H'o := {{dD)T) x H i{{dD)T) be the dual space of Hq with the duality between Hq and H'q 

defined by 

= {^1, V'l) + (V' 2 , ^^2) for (V?!, ip2f G H'o, {Yi, Y2f G Ho- 

In the following, we first prove the invertibility of A : Ho —>■ H'o- Define a bilinear form a(-, •) on Ho x Ho 
associated with the operator A by 


((ui, U2f, {vi, V2f') ■-= ) 


= (Viiui, vi) + ((i/ - iGii)u 2 , vi) + {v 2 , {-^I + Nii)ui) 
+{V 2 , W 11 U 2 ) + {V 2 , XU 2 ) for (mi, U 2 Y', {vi, V 2 )'^ G Ho- 
There exists a positive constant Ci such that 

|a ((ui, U2)^, (?^l, V2)'^) I < Ci\\{ui, U2 )^\\-Ho IK^li V 2 )'^\\-Ho- 


(2.4) 


Note that 


a ((mi, M2)^, (wi, U2)^) 

= {ViiUi, Ui) + {-K11U2, Ui) + {U2, +N11U1) + {U2, W11U2) + {U2, XU2) 


fVn -iLii 
liVn IDii 


where the operator 


Aq 


{U 2 , XU 2 ), 


Vn -Hn 
Hn Wn 


is positive on Ho (see ll|, Theorem 3.11]). Then we have the estimate 

a ((Ui, U2f, {Ui, U2f) > C2\\{ui, M2)^||«o, 


(2.5) 


where C2 is a positive constant. From (12.41) and (12.51) we conclude that A : Ho H'o is an isomorphism. 

In addition, it follows from Corollary 3.14 in ll| that 

h + K22-- H^’i{{dn)T) ^ H^'i{{dn)T) 

is also an isomorphism. Since r(a:, t; y, s) is smooth enough for Xy^y,0<s<t<T, the following operators 


K 21 ■- H^-^{{dn)T) H^^^{{dD)T), IT 21 : iL2>4((50)T) ^ H-2’-^{{dD)T), 

C 12 : H-i^-i{{dD)T) ^ Hhi{{dn)T), K 12 : Hhi{{dD)T) ^ (((90)t) 


are compact. Therefore, the system (12.21) is Fredholm with index zero. To show the unique solvability of 
(1^ . it suffices to show its uniqueness. 

Let {(fii, (p 2 , '-po)'^ G "H be the solution to the homogeneous form of (12.2p . That is. 


Dll 

kl-Kii 

K 21 \ 

1 Pi 

+ N^ii 

ITii + XI 

-fD21 

Pi 

D 12 

-Ki2 

+ K 22 J 

\P3 


O' 

'^2 I = I 0 


( 2 . 6 ) 
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Define 


j(x, t) := - f f r(x, t- y, s)(pi{y, s)d(T{y)ds + f f ^\ 2 {y, s)da{y)ds 

Jo JdD Jo JdD 9 ^\y) 


fl 

Jo Jdi 


5r(a;, t-, y, s) 


lo Jdn 9u{y) 
Note that w{x, t) satisfies 


ipoiy, s)da{y)ds, (x, t) G (M” \ {dD U 90))^ . 


dtw — Aw = 0 in (R" \ 0 )t 


and 


w = 0 at t = 0. 

Letting x tend to the boundary dfl, we derive from the jump relations of layer potentials that 


„+ — 


1 


= —Vi2(pi + Ki2(p2 — (t;! + K22)^3- 


(2.7) 

( 2 . 8 ) 

(2.9) 


Throughout this paper we use ‘+’ and ‘—’to denote the limits taken from the exterior and interior of a 
domain, respectively. It implies from the third equation of (12.61) that 


w~^ = 0 on {dil)T- (2-10) 

Then, using the same argument as in [i^, Appendix A], we can prove that w = 0 in (R" \ Q)t- It follows 
that dyW'^ = 0 on (cAI)^, and hence dyW~ = 0 on (i9D)t due to the jump relations of layer potentials. 

On the other hand, the function w defined by (1231) also satisfies the heat equation in (D \ D)t- Using 
the jump relations of layer potentials again, we derive that 

w = -Uiitpi + Kiiip 2 + - K2i<f3 on {dD)T, 

dvW = - Niiipi - Wii(p 2 + W 2 W 3 on [dD)T- 

In terms of the first two equations of (12.6L we have w = ip 2 and d,^w = \(p 2 on {dD)T- In conclusion, the 
function w defined by dlH) satisfies 


nSt - A)w = 0 in(D\D)T, 
I duW — Xw = 0 on {dD)T, 

I dvW = 0 on (9 D)t, 

[ w = 0 at t = 0. 


( 2 . 11 ) 


By the uniqueness of solutions to dm, we have w = 0 in (D \ D)t, and hence 

w~ =0 on (9 D)t. (2 -12) 

Combining (12.101) and (I2.12p . we obtain (/Ja = w~ — = 0 on (90)^. As a consequence, we have 

/ V,, ^ ( Q\ 

V-i/ + A^ii Wii + Xl) \y>2) Voy 

Now by the invertibility of A, we conclude that = (^2 = 0 on {dD)T- Due to the Fredholm theory, the 
unique solvability of the system (12.21) is justified. Moreover, according to the mapping properties of heat 
layer potentials, the function u expressed by dm is the desired solution in 77 ^’2 ((D \ D)t)- The proof is 
complete. □ 


Based on our Theorem 1 2.1 1 and Lemma 2.4 in 


ll| . we can define the Neumann-to-Dirichlet map Ad by 


Ad : 77-^-i((5D)t) ^ 775.I((5D)t), / ^ 


(2.13) 
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If there is not any cavity inside ft, i.e., D = %, the unique solvability of the forward problem in H^’i(ftT) 
can be found in [ll|. In this case, we denote the Neumann-to-Dirichlet map by A 0 . Taking the Neumann-to- 
Dirichlet map as the measured data, our inyerse problem for is formulated as follows: 

Inverse Problem: Reconstruct D from Ajj. 


3 The linear sampling method for the inverse problem 

In this section, we present a linear sampling-type method for the inverse problem formulated above and give 
its mathematical justification by using the layer potential argument. Let t) := r°(a;, t; y, s) be the 

Green function of the heat operator in fix with Neumann boundary condition on its boundary {dfl)T- Then 
the linear sampling method for the heat equation is based on the characterization of the approximate solution 
to the Neumann-to-Dirichlet map gap equation 

(AD-A^)g = T°^y^^^{x,t), {x,t)G{dn)T, (3.1) 

where s S (0, T) is a fixed time and y £ ft is the sampling point. By this characterization, we can define an 
indicator function as a mathematical testing machine to reconstruct the location and shape of D. 

To begin with, let us define the operators S, H, A and F as follows. 

• Define S' : ((SO)t) ^ and H : {{dft)T) H-h-i{{dD)T) by 

S : / !-)■ u^\dt, (3-2) 

and 

H : - Xu-^)\(^od)t: (3-3) 

respectively, where is the solution to 

( {dt — A)u^ =0 in ftx, 

< = f on {dft)T, (3.4) 

I = 0 at t = 0 . 


• Define A : FI 2 ^ i{{dD)T) ^ ft) t) hy 

^4 : 5 ( 3 . 5 ) 

where z® is the solution to _ 

({dt-A)z3=0 m{ft\D)T, 

I - Xz^ = g on (SD)t, 

I d^z9 =0 on {dft)T, ^ ’ 

[ z® = 0 at t = 0. 

• Define F : F[~^'~^ ((9D)t) —>■ Fl^’i ((i9II)t) by F := Ad — Ag. 

To characterize the solution to (EH), we first investigate the operator F and prove the following lemmas. 
Lemma 3.1 The operator F can he factorized as F = —AFl. 

Proof. For / G H~^~^{{dft)T), let u and v be such that 

f{dt-A)u = 0 in {ft\D)T, 

I dvU — Am = 0 on {dD)T, 

I d^u = f on {dft)T, 

[ u = 0 at t = 0, 
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and 


{dt — A)v = 0 in fi-r, 
diyV = f on {dn)T, 

V = 0 at t = 0, 


Then it holds that 


respectively. Define w as the solution to 

{dt — A)w = 0 in (D \ D)t, 
dvW — Xw = —{djjV — Xv) on {dD)T, 
dvW = 0 on {d^)T, 
w = 0 at t = 0. 

{dt — A){u — V — w) = 0 in (D \ £>)t, 
d^{u — V — w) — X{u — V — w) = 0 on {dD)T, 
dv{u — V — w) = Q on (9 D)t, 
u — V — w = 0 att = 0. 

The uniqueness result in Theorem 12.1 1 savs that u — v = w in {{XI \ D)t)- Hence, we have 

A{-Hf) = w;|(an)T = {u- v)|(an)T = - A^)f = Ff, 

which completes the proof. □ 

Lemma 3.2 The operator H : {{dXl)T) H~i’~i{{dD)T) is continuous and has a dense range. 

Proof. We prove the result using the layer potential argument. For / S {{dTl)T), we consider the 

following initial-boundary value problem: 


{dt — A)u = 0 in XIt, 
dvU = f on {dXl)T, 
u = 0 at t = 0. 


(3.7) 


Express its solution by a single-layer heat potential 


{x, t) = Votp := / r(a;, t; y, s)ilj{y, s)da{y)ds, {x, t) S XIt 

do dan 


(3.8) 


with an unknown density tp £ H i{{dXt)T). Using the jump relations of layer potentials, the problem 
(E3 is reformulated as the boundary integral equation 


{-I -f N 22 )iI^ = f on {dn)T, 


(3.9) 


where the operator i/ + N 22 '■ H 2 ’ 4 (( 90 ) 7 ') H 2 - 4,{{dXl)T) is an isomorphism (see Uj, Corollary 
3.14]). 

Define UV' := {di,{Vo'tp) — X{Votjj)) \{dD)T- Then the operator H can be represented as 

H = V{^I + N 22 )-\ 

Thus, the continuity of H is evident. To prove the denseness property of H, it suffices to show that the 
operator V : H~2’-4{{dXt)x) —>■ H~^'~4{{dD)T) has a dense range. 

Indeed, by direct calculations, we have 


Vil){x, t) = 


0 JdQ. 


M{x, t] y, s)tj}{y, s)da{y)ds, {x, t) £ {dD)T 


(3.10) 



with M{x, t; y, s) := 3j,(a;)r(a;, t] y, s) — A(a;)r(a;, t; y, s). Denote by V* : H^'i{{dD)T) {{dfl)T) 

the transpose of V in the sense that (D'i/'i v) = (V’; ^*v) for any ip G H~^~i ((90)t) and 77 G i?(( 9 £))y). 
It follows that ^ 

y*v{y^ s) = f f M{x, t; y, 5 ) 77 ( 2 ;, t)da{x)dt, (y, s) G (5 D)t. (3.11) 

J s J dD 

Then, to show that the range of V is dense, we are led to prove the injectivity of V*, that is, 77 = 0 in 
H^’i {{dD)T) if we have ¥*7] = 0 on (917)^. To this end, we define 


w{y, s)= [ [ M{x, t; y, s)y(x, t)da{x)dt, (y, s) G (M" \ dD) 

J s J dD 

Noticing that w satisfies 


(3.12) 


dgW + Aw = 0 in (R" \ O)^, 

Til = 0 on (9 D)t, 

Ti = 0 at 5 = T, 

we have tc = 0 in (R" \ D)^ by the same argument as in [iH Appendix A]. Using the unique continuation 
principle, we further have 


and therefore 


Note that 


dw+ 

dv 


= 0 in (R" \ D)t, 

(3.13) 

0 , w'^ = 0 on {dD)T- 

(3.14) 


w{y, s)= / M{x, t] y, s)ri{x, t)da{x)dt 

J s J dD 

= / {du{x)^{x, t; y, s) - A(a;)r(2;, t; y, s)) r]{x, t)da{x)dt 

J s J dD 

= / {duix)^iy, T - S] x,t)- A(a;)r(y, T - s; x, t)) 77 ( 2 ;, T)da{x)dr 

Jo JdD 

where 77 ( 2 ;, r) = 77 ( 2 ;, T — t). Using the jump relations of layer potentials, we have 

dw'^ dw~ 


dv dv 

It implies from (13.141) and (13.151) that 

dw 


dv 


= A 77 , Ti"*" — w =77 on {dD)T- 


— Xw =0 on {dD)T- 


(3.15) 


(3.16) 


Observe that w also meets 


dsW + Aw = 0 in Dt, 
w = 0 at s = T. 


Then, by the uniqueness of solutions to the backward problem, we obtain that 


w = 0 in Dt- (3-17) 

Thus, it can be concluded from (13.131) . (I3.15p and (13.171) that 77 = 0. This completes the proof. □ 

Lemma 3.3 The operator A : H~^~i{{dD)T) — t Hi i{{dV,)T) is injective, compact and has a dense range. 
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Proof. The injectivity can be easily seen from the unique continuation principle for the heat operator dt — A. 
We now show the denseness. Let Qj G H~i ~i {{dD)T) (j G N) be such that the linear hull of {gj} is dense 
in 24 ((3 £))t). Then it is enough to prove / = 0 if 

/ ipjfda{x)dt = 0 
J{dn)T 

for / G {{dil,)T) and all (fj := Ag^ = Here we recall that is the solution to (13.61) with 

9 = Si- 

Let V G ((O \ D)t) be the solution to the following backward problem 

|'(9t+A)z; = 0 in {fl\D)T, 

I dvV — Xv = 0 on {dD)T, 

I dvV = f on 
[ n = 0 at t = T, 


and set Zj = . Then, we have 


0 = / (vAzj — ZjAv)dxdt 
A^\D)t 

/ {d,yZj V — duV Zj)da{x)dt — / {d^zj v — d^v Zj)d(j[x)dt 
J{an)T J{dD)T 

/ gjvda{x)dt. 
hdD)i 


So = 0 on {dD)T- By the boundary condition of u, we further have dvv\(^QD)^ = 0. Therefore, r; = 0 in 
(H \ D)t, and then / = d^v\(^oQ)^ = 0. 

Finally, let us prove the compactness of A. Indeed, there exists a unique density (/? G {{dD)T) 

such that the solution z® to (1^ with g € H 2 ’ 4 [[dD)T) is given by 

z9(x, t)= [ f r°(x, t] y, s)(p{y, s)da{y)ds. 

JQ JdD 

Notice that r°(a;, t; y, s) is smooth enough for x G Oft, y G dD and 0 < s < t < T. It can be concluded that 
Ag = z®|(an)T £ C°°{dV, x [0, T]), and hence A is compact. This completes the proof. □ 


We are now in a position to state our main results, which motivate the linear sampling method for 
reconstructing D. 

Theorem 3.4 Let s G (0, T) he fixed. For y € D, there exists a function g^ G {{dXl)T) satisfying 


such that 


and 


^{v,s)\\HFi((dn)T) ^ ^ 


^hm^||55^||^,4(^^)=oo. 


Proof. According to Lemma 13.21 for any £ > 0 there exists a function g^ € H 2 


Hgy + (- 






< 


||A|| 


(3.18) 

(3.19) 

(3.20) 


4 ((9H)t) such that 

(3.21) 
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where ||A|| is the norm of the operator A defined by (13.51) . Since satisfies the heat equation in (f2\ D)t 

and = 0, it holds that 




By Lemma ixn we have 




< 






11 £. 
H^2’ 4((dD)T) 


Hence, due to the boundedness of iF, we have 


Hr 




> ciiiFri 


H^^’^S{(dD)T) 


> c 

> c 


9<2^ly,s)->^^ly,s) 




Hs’ + ( 








Ce 




(3.22) 


where C is a positive constant. 
We next prove 




11 —)► oo as y ^ dD. 

H^^-S{{dD)T) 


(3.23) 


Since r°(x, t; y, s) — r(a;, t; y, s) G C°°{flT), that is, r° has the same singularity as F at (x, t) = {y, s), we 
only need to prove that 


\d,T 


iy^s) 


-XT 


(y, s) 


iH-h-iaoDM^^ ^y^dD. 


(3.24) 


We consider here the case for n = 3. Under the assumption that dD is of class C^, for any point xq G dD 
there exists a C^-function $ such that 


D (~l B{xo, r) = {x G B{xo, r) : X 3 > d>(a;i, X 2 )} , 

where B(xo, r) is the ball with radius r > 0 and center at xq. We choose a new orthonormal basis {ej}, j = 
1, 2, 3, centered at xq with 63 = —v, where v is the unit outward normal vector to the boundary at xq. The 
vectors ei and 62 lie in the tangent plane to dD at xq. Let x be the local coordinate defined by the basis 
{cj}. We introduce the local transformation of coordinate rj = F{x) as follows: 

r]' = x', T]3 = X3- $(x') with x' = (xi, X2), rj’ = (? 7 i, m)- 

We note that <i>(0) = Va;'4>(0) = 0. For x G dD and y G D, let rj := F{x), ^ := F{y) with y = ( 77 ', 0), 
^ = (0, ^ 3 ). Then it holds that 

|x - y\^ = \F-\y) - F-\0? < C 2 |r? - = 02 (^ 3 " + bf )• 


We now pick up the dominant part of ()3.24|) . that is, 


l|5i?3r(«,i.) - AF(j_^)|| i,_i 


H 2’ 4((D2)t)’ 


(3.25) 


and ignore all the other terms which are bounded as ?/ —>■ dD, where D 2 ■= [—1, 1] x [—/, 1] with 0 < Z <C 1. 
To estimate (j3.25L without loss of generality, we assume s = 0 and introduce an auxiliary function 


t) 


cxW. t)t “e 
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with 0 < a < 1/2, where xiv 'j ^ smooth cut-off function such that it vanishes near t = T and = 
±l{j = 1, 2) but equals one in a neighborhood of the origin (ry', t) = (0, 0). Let us first show that G 
((D 2 )t) C H2’4 ((D 2 )t) and n n ^ sufficiently small constant c. 

Indeed, by direct calculations, we have 


‘Pt =cx{v',t)e (^-at “ ^ 

Pv, = -^cc2x(?y', 


lr?'P \ I 'G 

+ C2^t-“-M +cxt(ry', — 

Vj + cXvjW^ 


So, for 0 < t <C 1, and can be estimated by 


Note that 


Then we have 


10 JD 2 

In addition, we deduce that 


\pt\<C3t “ , \iPr,j\<C3t “ 


dr]'= [ e-^^^dr]' = 0{t). 

[ [ \^ 7 ]j\‘^dxdt < C 4 j tdt = C 4 f t~‘^^dt<oo. 

Jo JD 2 Jo Jo 




1/2 


< 


Vr(3/4) 

y (t - ^ \dTp{x, T)\^dx^ dr 


r(3/4) 


< 


C5 / {t — x) 
Jo 




(3.26) 


(3.27) 


which guarantees the integrability of in under our condition 0 < a < 1/2. Thus, we conclude that 

ip G H^’4 ((D 2 )t) C H2’4 ((D 2 )t) and i i,,„ , < 1 by taking a sufficiently small constant c. 

/f2’? {{D 2 )t) 

Next, we compute the norm of driX{i,s) ~ g) using the duality. In fact, we have 

||^i?3r({,s) - 11^- 1 ,_ 1 


> 


C7 


> 


cs 

87r3/2 

Cs 

87r3/2 


y y 3/2e ^cxt “e ^ dry'dt 


-l/2-a -C2 




■ ^ 2 ) 

/«i T“"5e-'==idT 00 as ^ 3 ^Q {y ^ dD). 


Now we can conclude from (13.221) and (13.231) that the blow-up properties (13.191) and (I3.20p hold, 
proof is complete. 


The 

□ 


To further investigate the behavior of the density as y approaches to dD from the exterior of H, we 
need the following lemma. 

Lemma 3.5 For any fixed s G (0, T), the Green function T®^ s)(®’ range of A if y G fl \ D. 
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Proof. Suppose that s)(a^! t) is in the range of A. Then there exists a function g € H i’ i{{dD)T) 
such that 

i«®l(cin)T = r°y_s)l(an)T> (3.28) 

where w® G ((O \ D)t) is the solution to 

{dt — A)'w^ =0 in (f] \ D)t, 
dvW^ — Aw® = g on {dD)T, 
duW^ = 0 on {dVL)T, 
w® = 0 at t = 0. 

From (13.281) and 9,.w®|(an)T = = 0; it follows that w® = in (fl \ (£> U {j/})).^, and hence 

^ ^ 2 / G \ £> can be either y e dD ot y £ fl\D, and for 


each case it holds ||r°^ 




= oo. This gives a contradiction and completes the proof. 


□ 


In contrast to Theorem 13.41 for y £ we now establish the following blowup property of the density 5 ® 
for y ^ D. 

Theorem 3.6 Fix s £ (0, T) and let y £ D, \ D. Then, for every e > 0 and (5 > 0, there exists a function 
3l,s ^ H-^-^{{dTl)T) satisfying 


such that 


and 


\\^9e,s ^(y.s)ll_H-i'3((an)T) ^ 


lim llg® oil 11 


lim ||S'g® aII -1 1 , 


= 00. 


(3.29) 

(3.30) 

(3.31) 


Proof. By Lemma [3.31 for arbitrary i5 > 0 there exists a function fg £ H 2 -- 4((9£))t) such that 

Since F^^ is not in the range of A due to Lemma 13.51 we have 

as <5^0. 


(3.32) 


Recall that the range of Ft is dense in iL i - 4 (((5£))t)- We can find a function g^ g £ H 2 ’ 4 ((9n)T) such 
that 

nV j- fy\\ , , ^rl(\\A\\J-^\ ( 3 . 33 ) 




and hence 

So, it can be derived that 

r,V vO 


Ufs+AHglJ.. 


e.»"_H-2’4((an)r) 


< e. 




H^’t{{an)T) ^^9e,s i^(!/.s)lli/i.3((an)T) 


< £ + (5. 

Moreover, we obtain from (13.321) and (13.331) that 

\\^9ls\\H-h^liiaDM^^ asd^O, 
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which implies 


Il5" 




—>■ OO, 


IIV 


e,i5l 




OO as (5 —)► 0. 


This completes the proof. 


□ 


By Theorems 13.41 and 13.61 Sg^ can be taken as an indicator function for reconstructing the boundary of 
the cavity D, where is the approximate solution to (13.11) for the sampling point y G fl. 


4 The asymptotic behavior of the indicator function 

In this section, we show a short time asymptotic behavior of the indicator function Sg^ as y approaches to the 
boundary dD from the interior of D. To do this, we analyze the asymptotic behavior of the reflected solution 
of the fundamental solution. As we know, the reflected solution is the compensating term of the Green 
function for the related initial-boundary value problem. This means that we actually provide a pointwise 
short time asymptotic behavior of the Green function near dD. From our asymptotic analysis, we can 
establish a pointwise reconstruction scheme for the location and shape of the cavity, and can also know the 
distance to the unknown cavity when we probe it from the interior of D. Without loss of generality, we only 
consider the three dimensional case. 

Let t) be the reflected solution of t). That is, 

( {dt - A)r(y_s)(x, t) = 0 in Dt, 

- A)r(y_^)(x, t) =t) on {dD)T, 

I t) = 0 ioi X G D, t < s. 

Then, by the well-posedness of the above initial-boundary value problem, we obtain from (I3.21|) that 

< C'le, (4.1) 

where Ci is a positive constant. Denote by s)(x, t) the reflected solution of r(y s)(x, t), i.e., 

( {dt - s)(x, t) =0 in Dt, 

< {d^ - \)R(^y^s){x, t) =-{d,y-X)T(^y^^){x, t) on {dD)T, (4.2) 

[ i?(y^s)(x, t) = 0 ioT X G D, t < s. 

Since T^^ ~ r( 2 /.s)( 3 ^) t) for {y, s) G Dt is smooth in Dt, there exists a positive constant C 2 such that 

ll^(y, s){x, t) — s)(x, t)||^ ^ 2 - (4-3) 

Combining (14.11) and (14.31) yields 

ll(S'g^) (x, t) - R{y,s){x, < Cie + C 2 . (4.4) 

This implies that Sg^ and i?(y g) have the same blow-up behavior as y approaches to dD. Hence, to get the 
asymptotic behavior of Sg^, we only need to show the corresponding result of R(y^s){x, t). 

To begin with, we locally flatten the boundary dD. For any point z G dD, there is a diffeomorphism 
^ which transforms z to the origin 0 such that 

<h(D) C R3 := ^ 2 , 6) e R' : 6 < 0} . 

Without loss of generality, we assume 0 G dD with z/(0) = 63 and locally flatten the boundary dD around 
0. Let ^ = 4)(x) and 77 = $(?/). Denote by J(x) the Jacobian matrix of <1>. We have 

J(x) = V<l>(x), J(0) = I, i^ix) = 
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Set R{^, t; r], s) := R{x, t\ y, s) and r(^, t; y, s) := r(a;, t; y, s). We deduce from (14.21) that i?(^, t; rj, s) 
near 0 satisfies 

r {dt - V5 • M{OV^)R{t t- r;, s) = 0, t€ (0, T), ^3 < 0, 

1 63 • t; 77 , s) - v, s) 

S = - (e3 • (M(e)V5f)(e, t; 77 , s) - HO\m t', V, s)) , 

I R{^, t] 77 , s) = 0 , t< S, ^3 < 0, 
where M(^) = (Let t; 77 , s) be such that 


t € ( 0 , T), ^3 = 0 , 


( 4 . 5 ) 


{dt - A«)f^+(e, t; 77, s) = 0, < e (0, T), ^3 < 0, 

(% - A( 0 )) W+(e, t- 77, s) = - (% - A( 0 )) f (^, t; 77, s), t e ( 0 , T), ^ = 0 , ( 4 . 6 ) 

t; 77, s) = 0, t < s, ^3 < 0. 


In the following, we first estimate the difference R{^, t; 77, s) — t\ 77, s), and then analyze the 

asymptotic behavior of W+(^, t; 77, s). As a consequence, the asymptotic behavior of R{^, t\ 77, s) is obtained. 


Lemma 4.1 Let 77 = (771, 772, 773)^, where 773 = —e /2 aTid I77I < ce for small e > 0 and some constant c > 0 . 
Then there exists a positive constant C such that 


R{y, s + e^; 77, s) - W+{y, s + e^; 77, s) 


< Ce-2 


as e —>• 0. 


( 4 . 7 ) 


Proof. Define P{f, t; 77 , s) := R{f, t; 77 , s) — W+(^, t; 77 , s). We derive from (j4.5p and (14.61) that near 0 the 
function P(^, t; 77 , s) satisfies 

(a* - Aj)P(^, t; 77 , s) = V • (A/ - /)Vi?, t G (0, T), 6 < 0, 

(aj 3 - A(0)) P(e, t- 77 , s) = (a^ 3 (i? + f) - 63 • [A/(0V^(I? + f)]] 

-(A(0)-A($-i(e))|J^Kai) (^ + f), t G (0, T), 6 = 0, 

-P(^, s) = 0 , t < s, ^3 < 0. 


Denote by f(,,,s)(C; t) '■= f(C: t] s) the solution to the following problem; 

f -{dt + A^)f(^, t- 77 , s) = d{t - s)S{^ - 77 ), t G (0, T), ^3 < 0, 

I (% - A(0)) f (C, t; 77 , s) = 0, t G (0, T), ^3 = 0, (4.9) 

[ f (C, V: s) =0, t> s, fs < 0. 

Let Q be a small bounded domain in the lower half-plane such that 0 G dQ (~l SR?. C d{^{D)) and Q totally 
lies in $(D). We decompose the boundary of Q into two disjoint parts such that dQ = dQi U dQ 2 with 
aQi = dQ (~l aK^. Since we only consider the singularity locally at 0 , we have 


Pin, s) it t) = - 


L Lq 'r)a^(^)f(5_t)(2;, t) - d^(^,)P(^r,,s){z, r)f(5_t)(2;, r)) dcr{z) 

^ Is I dQ ~ ^ 

~ I “ ^) s) {z, t) ■ vf t) {z,T)'j dz 

- dr (P(^r,,s){z, T)a^ 3 f( 5 ,t)( 2 :, t) - d^^Pin, s){z, T)t(^^^t){z, t)) dz’ 

Js JdQi ^ ' 

+ [ dr f (e 3 -{M{z) - I)X/R(^r,,s){z, t) ■t(^^^t){z, t)) ds{z) 

Js JdQi ^ ' 

- Is dr j (^{M{z) - I)X/R(^r,,s){z, t) ■X/r(^^^t){z, T)^dz + 0{1). 


( 4 . 10 ) 
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where fj = (771, 772, and Cj (j = 1 , • • • , 4 ) are positive constants. 

On the plane ^3 = 0 , it holds that 


P{n,s)iz, T)dz3t{i,t){z, t) - dz^P(n,s){z, T)f(5_t)(z, t) 

(a,3(i? + f) -63 • [M(z)V,(i? + f)]) r) 

^ (A( 0 ) - A($-i(z))| J^K^)I) (R + r) f (5, t)(2, r) 

(e3-[(/-M(z))V,(/? + f)])f(5,,)(z, r) 

[A( 0)(1 - \J^Yz)\) + (A( 0 ) - A($-i(z)))|J^K^)|] (i? + f)f(5,,)(^, r) 


+ 


+ 


< Co\z'\(^VR^r,,s)iz, t) + Vf(^_3)(2;, r) + Yr,,s){z,T) + t^r,,s){z, t)^ f(5_t)(z, t) . ( 4 . 11 ) 

Set z' := (zi, Z2), rj' := (771, 772) and C '■= (Cii ?2)- Using Lemma 3 in [l 7 |, Chapter 1 ], we derive that 


< 


f dr f \z'\( WR(^^^s){z, t) + VT(^^^s){z, t)) f(5_t)(2, t) 

Js JdQi ^ ^ 

Cl f dr [ {\z' - 'n'\ + \'q'\) (wYr,,s){z,T) + wr(r,,s){z, t)) t(^^^t){z,T) 

J s J '' ' 


< 


IdQi 

Xo. I"'-”'1 (-s)'’) 


dz' 

-eT+4" 


,6. /“./■, 7,. ^-2 ^ \z'-i?+vl\u ^-3/2 r k' 

+C'3 J dr \v\[r- s) exp -— j {t - t) ' exp - 


< 


+ C'5|»7'I%^^3’^^ 


/OQi 


-?f+g 

C2(t - r) 
C7(t - r) 

C2[t - t) 

-s)J ' 


dz' 


( 4 . 12 ) 


Noticing that 


f dr [ \z'\( R(^rj,s){z,T) + r(^r,,s){z,T)) f( 5 ,t)( 2 ;, t) 

Js JdQi ^ ' 


dz'. 


has a less singularity than ( 14 . 121 ) . we omit the estimate of this integral here. 
In addition, the integral 
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can be estimated in the same way as in (14.121) . As for the volume integral in (I4.10p . we have that 


dr {[M(z) - I)VR(^ri, s) ( 2 , r) • Vf () (z, r)) dz 

dz 


.t . 



< Cg / dr |z| 

VP(^,s)(z, r) 

Vf( 5 ,t)(z, r) 

Js Jq 


<Cg{t-s) ^exp(—+ CiolryK^ - s) 


C 4 (t — s) 


<C9{t-s) ^exp ()+C'io|? 7 |(t-s) 3/2gxp(- 


C 4 (t — s) 

16 +»?3P 


C 4 (t — s) J 


(4.13) 


Take ^ = rj with I 77 I < ce, 6 + ^3 = ~£i t — s = e^. Then we finally get the following estimates for e —>■ 0: 


f dr f T)a3f(j_t)(z, r) - a3P(^_^)(z, T)f(^_t)(z, t)) dz' 

Js JdQi ^ ^ 


< Ciie- 


z) - /)VP(^_^)(z, r) • Vr(^_ 4 )(z, t) \ dz 


< Ci2e 


which lead to the estimate (IQ) by (|4.10l) . The proof of this lemma is complete. 


□ 


Next we derive the expression of t; rj, s), and then show its asymptotic behavior. 

Lemma 4.2 The solution t; rj, s) to (14.611 can be expressed by 

VL+(6 t- r,, s) = e^^dT [ e^^'<' f + ^ 

IGtT'^Z Jcr-ioo JR2 0(0 - Ao) 

where 6 := (6j 6)1 v' •= (dii V2), Aq := A(0) and 0 := a/t + |C'p with Re0 > 0. 

Proof. Define id(6 t', Vj s) := lT+(5, t; rj, s) + f(^, t; rj, s). We know from (14.61) that id(6 t', s) satisfies 


{dt - A5)id(5, t; 77 , s) = S{t - s)S{^ -rj), t€ (0, T), ^3 < 0, 
(% - A(0)) P(5, t- v,s) = 0, tG (0, T), ^3 = 0, 
t', V, s)=0, t< s, ^3 < 0. 


(4.15) 


Denote by H the Laplace transform of H with respect to t. Then we have 

f (r - Aj)id( 5 , t- r], s) = e“^®(5(^ - p), ^3 < 0 , 

j (% - Ao)id(6 d, s) = 0, 6 = 0. 

We now look for the solution to (I4.16|) in the form of 

^(6 d, s) = id =^(6 t] r], s), ±(^3 - 773 ) > 0 , 

where .ff^(6 d: s) satisfy 

{ {t - A^)id±(^, t; rj, s) = 0 , ±(^3 - 773 ) > 0 , 

\ H+-H- = 0, ^3 = V3, 

) %(P+ - H-) = -e—<5(6 - 6), 6 = d3, 

[ (% - Ao)id+(6 6 d, s) = 0, 6 = 0. 


(4.16) 


(4.17) 
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Denote by the Fourier transform of with respect to ^' = (^i, ^ 2 ) and let = (Ci, C 2 ) be the Fourier 
variable corresponding to . Then we have 


(t+ IC'P - i9|3)</5=^ = 0, ±(^3-7?3)>0, 

ip+ -ip~ = 0 , ^3 = 773 , 


,+ _ ,n-\ = -p-rSp-iC'-v' 


) = -e 


^ ?3 = Vs, 


(4.18) 


(% - Ao)(/3+ = 0 , ^3 = 0 . 

Since the operator — 0^ has the fundamental solutions we seek the solution to (14.181) in the form 

of 

ip'^ = + C2 , ip~=ce^^^. 

The transmission conditions at ^3 = 773 give 


0 = ip~^ — ip = cf + c^e ®’'^ — ce®’’^, 
_e-(^*+<'-V) = 0c+e®''3 - Oc^e"®''" - 0ce®''T 

The boundary condition at ^3 = 0 leads to 

0 = (0C+ - 0C+) - Ao(c+ + c+) = (0 - Ao)c+ - (0 + Ao)c+ 

which implies 


(4.19) 


0-A, 


C 9 — 


0 + 

-c7. 


0 + Ao ' 


(4.20) 


From (14.191) and (14.201) . we obtain 


q -20 e , C 2 -^0 e e 


Thus, ip'^ is expressed by 




26- r e-A, ■- J- ('‘-^1) 

We note that in the brace of (14.211) the second term comes from the fundamental solution, while the first term 
corresponds to the reflected solution VF+. So (14.141) is obtained by taking the inverse Fourier and Laplace 
transforms of (14.211) . The proof is complete. □ 

Finally, we show the pointwise asymptotic behavior of t\ 77 , s). 

Lemma 4.3 Let rj = ( 771 , 772 , 773 )^, where 773 = —e/2 for small e > 0. Then we have 

— 3 —3 

W+( 77 , s + e^; 77 , s) = ~ 8 ^ 174 ^ 3/2 ^ as e ^ 0 . 

Proof. To analyze the asymptotic behavior of W~^ , let us first calculate 

/•a+ioo ft I \ 

/ = /(C', t, s) := / 6, m < 0 

J <7—ioo vJ(^0 Aoj 


(4.22) 


(4.23) 


by converting the line integral from a — ioo to tr + zoo into a closed contour so that we can apply the residue 
theorem. The contour used here is described in Figure HTTl 

The contributions from the arcs ABC and FGH are negligibly small as i? —?► 00 by Jordan’s lemma. Let 
Fp be the circle with radius p centered at (—|CTj 0)- Then, by setting r + = pe*®, we can show that 

[ p(t-s)T Q + Aq e(;3+?73)4- 

^ 0(0 - Ao) 

— f _ v^e / + Aq _ dO ^ 0 as o —>■ 0 (4 241 
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Figure 4.1: In the contour marked with arrow, R is the radius of the half circle and Fp is the circle with 
radius p centered at (— 0 )- 


„-(t-s)(’-+IC'P)_^ 

fl-icr iy/riiy/^ - \o) 

rfl-ICT 

f e 

0 






ir as i? —>■ oo. 


(4.25) 


Let R be the radius of the half circle in the contour. By taking r = —— t , we have 

[ At-s)T Q + Aq ei^3+ri3)p~ 

0(0 - Ao) 

-(t-s)(r+\c\^) iVr + \o ^iVF(; 3 +rj 3 )( 

i\/r(j\/r - Ao) 

^-(t-a)r -r + 2iXo^ + Ag 

Vr{r + Xl) 

2 iAoVF + Ag 

Jo V^ir + Xl) 

Similar to the derivation of (I4.25F we also have 

f Jt-s)r Q + Aq e(;3+r;3)p>, 

JeP 0(0 - Ao) 

^ r e-it-s)r as R 

Jo V^ir + Xl) 

Notice that (Ag — IC^Pi 0) is contained in the interior of the contour. By the residue theorem, we get 

-2ze-(‘-«)l«'I^Re ( T . (4.27) 

\Jo V^v + \ ) / 

Define 

L = L{t, s, 6, m; Ao) := Re ( H . 

\Jo Vr{r + Ag) J 

Then W'^ can be represented by 


—?► oo. 


(4.26) 




— f p-^n(*~'S)+Ao(g3+?;3) _ _!_r 

27r2 SttS 




(4.28) 
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By direct calculations, we obtain that 


/R Jr 

= {y/t - [ e~^^^'>'‘dCi 

Jr 


Similarly, we have 


Jr 

Therefore, it follows from f|4.28p that 


In the sequel, we compute L. By setting p = y/r, we have 
u,,,2-p^ + 2zXop + Xl 


L = Re 

= Re 


p(p2 + Ag) 


^e*P(«^+''^)2pdp 


r+oo 2 I \2 

-p + Ag 


(4.29) 


p2 + A2 


/ +00 n 

, (4.30) 

-oo P +Ao 


To compute the integrals in (I4.30p . we introduce the contour described in Figure 021 so that the residue 
theorem can be applied. 



Figure 4.2: The contour marked with arrow contains the point (0, —Ag) in its interior. 


Note that 


-{t - s)p^ + ip (^3 + ? 73 ) = - ( Vt-sp - + r? 3 )^ - 


and define 


p := y/t-sp - ^ J — (6 + %)• 


2y/t^. 
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It holds that 

-/ = -iii - + %) |(i - s)-^/^p + ^(t - s)"^(C3 + ?73)| - “ s)”^ 

= -i{t - s )"^/^(6 + m)p + ^(6 + - s)"^ - - s)"^ 

= A + ibp, 


where 


^■=\(^3 + mf{t-s) ^-p^{t-s) \ b:=-{t-s) ^/^( 6 +»? 3 )- 


Using the residue theorem and noticing that the integrals on Rep = ±i? go to zero as i? —^ oo, we obtain 
from (j4.30|) that 


L = + {t - 


2 r+°° 




(Ss+ns) 

—Aoe '‘U-b 


2 /* + 00 


~2 (A - XDh + 26p^(i - s) 


— oo 
-1 


(-A + A2)2 + 62p: 


J_^ i-A+xir+b^p-’ ^■ 


(4.31) 


Now take ^ = rj with ^3 = 773 = —e/2 and t — s = e^. Then b = e ^ and A = e — p^). By direct 
calculations, we have 


_A2+A4-6V + 


2A2p2 + 2e2A4 - A2/2 


{-A + Aq) 2 + 62p2 
{A — Ao )6 + 2bp^{t — s)“4 


p 4 + (1/2 + 2e2A2)p2 + (e2A2 _ 1/4)2 ^ 
- AqE^ + 1/4 


{-A + Ag)^ + 6 ^p^ 

Hence, it can be easily seen that 

.2 -A^ + X^-b^p"^ 


I 

J —( 


pi + (1/2 + 2e2A2)p2 + (£2^2 _ 1/4)2 • 

dp = -y/n + 0{e^), 


+ 00 


{-A + A2)2 + 62p2 
-1 


p 2 {A - Xl)b + 2bp^{t - s) 


-dp = 0 ( 1 ) 


, {-A + Ag)^ + b'^p^ 

for e —>■ 0. Consequently, we have 

L = 47rAo — + 0(1) as e —0. 

Thus, in terms of (j4.29|) . we finally obtain 

g-3 _ g-3 

lU+(7?, s + £2; ,y, s) = —rj^ + 0(e 2) ^ ^00 as e ^ 0. 

This completes the proof. 


□ 


Based on the above lemmas, we conclude that 1U+ is the dominant part of i? as e —>• 0. More explicitly, 
the pointwise asymptotic behavior of i?(^, t; 77 , s) can be stated as follows. 


Theorem 4.4 Let rj = (771, 772, 773)^, where 773 = —e /2 for small e > 0 . Then we have 

^-3 


R{r], s + e ; 77 , s) = 


+0{e-^) 


8gl/47r3/2 


8gl/47r3/2 


00 as e —>• 0 . 


(4.32) 


Remark 4.5 For any fixed discrepancy in (13.181) . combining (j4.4p with (14.321) . we finally have the following 
asymptotic behavior 

^-3 

(4.33) 


{Sgy){y,s + e) = -^-^j^ + 0{e ) as e ^ 0. 
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5 Concluding remarks 


This paper investigated the inverse problem of identifying an unknown Robin-type cavity inside a heat con¬ 
ductor from boundary measurements. The so-called linear sampling method was established to reconstruct 
the shape and location of the cavity. Based on our well-posedness analysis of the corresponding forward prob¬ 
lem, we gave rigorously a theoretical justification of this reconstruction scheme by using the layer potential 
argument. Further, we proved a short time asymptotic behavior of the reflected solution of the fundamental 
solution, and hence the asymptotic behavior of the indicator function used in this method was obtained. 
Since the reflected solution is the compensating term of the Green function for the related initial-boundary 
value problem, we actually provided a short time asymptotic behavior of the Green function. The asymptotic 
behavior naturally yields a pointwise reconstruction scheme for the boundary of the cavity. We would also 
like to emphasize that from the asymptotic behavior we can know the distance to the unknown cavity as we 
probe it from its inside. However, to establish a pointwise reconstruction formula for the Robin coefficient A, 
we need to carefully examine the lower order term where the information about the Robin coefficient should 
be involved. This is one of our future works. Also, we intend to utilize this asymptotic behavior to generate 
a good numerical performance of this reconstruction scheme. 
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